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Abstract. Dust particles observed in extrasolar planetary discs originate from undetectable km-sized bodies but 
this valuable information remains uninteresting if the theoretical link between grains and planetesimals is not 
properly known. We outline in this paper a numerical approach we developed in order to address this issue for 
the case of dust producing coUisional cascades. The model is based on a particle-in-a-box method. We follow 
the size distribution of particles over eight orders of magnitude in radius taking into account fragmentation and 
cratering according to different prescriptions. A very particular attention is paid to the smallest particles, close to 
the radiation pressure induced cut-off size Rpr, which are placed on highly eccentric orbits by the stellar radiation 
pressure. We applied our model to the case of the inner (< 10 AU) (3 Pictoris disc, in order to quantitatively derive 
the population of progenitors needed to produce the small amount of dust observed in this region (~ 10^^ g) . 
Our simulations show that the coUisional cascade from kilometre-sized bodies to grains significantly departs 
from the classical dN oc R~^'^dR power law: the smallest particles {R ~ Rpr) are strongly depleted while an 
overabundance of grains with size ~ 2Rpr and a drop of grains with size ~ 100-Rpr develop regardless of disc's 
dynamical excitation, Rpr and initial surface density. However, the global dust to planetesimal mass ratio remains 
close to its dN cx R~'^'^dR value. Our rigorous approach thus confirms the depletion in mass in the inner /? Pictoris 
disc initially inferred from questionable assumptions. We show moreover that collisions are a sufficient source of 
dust in the inner /3 Pictoris disc. They are actually unavoidable even when considering the alternative scenario 
of dust production by slow evaporation of km-sized bodies. We obtain an upper limit of ~ 0.1 for the total 
disc mass below 10 AU. This upper limit is not consistent with the independent mass estimate (at least 15 Mq) in 
the frame of the Falling Evaporating Bodies (FEB) scenario explaining the observed transient features activity. 
Furthermore, we show that the mass required to sustain the FEB activity implies a so important mass loss that 
the phenomena should naturally end in less than 1 Myr, namely in less than one twentieth the age of the star 
(at least 2 10^ years). In conclusion, these results might help converge towards a coherent picture of the inner 
(3 Pictoris system: a low-mass disc of coUisional debris leftover after the possible formation of planetary embryos, 
a result which would be coherent with the estimated age of the system. 
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1. Introduction 

1.1. The 13 Pictoris system 

The dusty and gaseous (3 Pictoris disc has been inten- 
sively st udied since the fi r st res olved image was obtained 
in 1984 ijSmith fc TerrileL Il984|) . This system is particu- 
larly interesting since it is still one of the best examples of 
a possible young extrasolar planetary system. It should be 
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stressed that consideri ng the estimated age of the systern , 
i.e. at least 2 10 ''years (jBarrado v Navascues et aLlll999|) . 
/3 Pictoris is no longer in its earliest formation stage and 
that if planetary accretion had to occur then it should 
already be finished. 

The (3 Pictoris disc has been observed in various wave- 
leng ths and has a radi al extent of at least 1500 AU 
(^Larwood fc Kalasl l200l)) . Due to obvious observational 
constraints (10 AU at (3 Pic's distance correspond to 
~ 0.5"), it is mostly the outer part of the disc beyond 
~30 AU that has been extensively mapped and studied 
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with t he highest spatial r esolution (the reader might re- 
fer to lArtvmowicj l)l997f) for a complete review.) The 
global picture is that of a relative central gap in the dust 
density, followed by a density peak around 100 AU and 
a slow decrease outwards. The exact density profile re- 
mains relatively uncertain since it strongly depends on 
both the assumed size distribution for the dust popula- 
tion and the grain optical properties. As a consequence 
the total mass of dust is also relativel y uncertain, but 
might be of the order o f 0.3 to 0.5 ((Li fc GreenbereL 
ll998HArtvmowiczL[l997|) . Note that "dust" mass estimates 
always strongly depend on the upper size limit consid- 
ered, and that observations do not constraint very well 
the abundance of the bigger dust particles where most of 
the mass is supposed to be. Observations show mostly the 
signature of micron-sized grains visible through scatt ered 
hght fe.g. lKalas fc Jewittl.ll995HMouillet et al.Lll997(l and 
thermal emission l|Lagage^rPantinl~ll994|) . Millimetre- 
sized gra ins are detected by millimetre wavelength pho- 
tometry ijGhini et all 1 199 it Iziickerman fc BeckhnI Il99,ll) 
and by resolved imaging in the submil limeter d omain but 
with a poor angular resolution (jHoUand et al.l [19981 . 

Actually no direct information is available for all ob- 
jects bigger than a few millimetres. Analytical estimates 
have shown that the observed dust cannot be primordial 
since its expected lifetime imposed by the rate of destruc- 
tive mutual collision s is much sho r ter than the estimated 
age o f the system l|ArtvmowiczL Il997t iLagrange et al.l 
Thus dust must be constantly produced within 
the disc. Candidates for producing this dust are suppos- 
ingly kilometer-sized planetesimals whi ch may generate 
dust e ither by ev a porat ion of volatiles l|Li fc GreenbereL 
[1998"; Lecavclier, '1998") or /and by coUision al erosion 
(■Artymowicz...l997;.Mouillet et al.Ul997; .Augereau et al.l 



200 1|) . In either case, the total m ass of the disc must be 
dominated by these parent bodies. lArtvmowicj l)l997l) es- 
timated t hat if a steady state coUisional law in dN cx 
R-^-^dR l|Dohnanvi Il969(l holds from the smallest dust 
grains to the biggest planetesimals, then one might expect 
at least 140 of kilometre-sized objects. But, as noted 
by the author himself, such extrapolations remain very 
uncertain. 

Scattered light observations have also revealed several 
more or less marked asymmetries in the outer disc (see 
iKalas fc Je-^ Il995'. for a detailed presentation). Some 
of these asymmetries arc believed to be due to the pres- 
ence of an embedded planet. It is in particular the case 
of the slight warp in inclination (~ 3°) of the disc's mid 
plane observed up to 80 AU. This warp has been suc- 
cessfully interpreted by the dynamical response of a plan- 
etesimal disc to the pull of a Jupiter-like object located 
at about 10 AU from th e star on a slightly inclined or- 
bit ijMouillet et al.[ll997t) . To extend this scenario to the 
dust disc moreover allows to repro duce large-sca le vertical 
asymmetries up to about 500 AU (j Augereau et al.. 2001) . 
The planetary hypothesis is reinforced by new asymme- 
tries evidenc ed at mid-IR wavelengths in the inner disc 
( Wahhai et al. ,200a) . 



1.2. The inner disc 

Nevertheless, these indirect effects of an hypothetic planet 
are detected much further away from the star than the 
planet's actual location. As indicated in the previous sec- 
tion, there is a strong lack of data for the region within 
10 AU which is probably the most interesting area in terms 
of presence of planets and planet formation. Most of the 
informations on this region has been indirectly inferred by 
fitting the Spectral Energy Distribution (SED) in the near 
and middle infrared. There seems to be a general agree- 
ment on the fact that the inner part of the disc is signif- 
icantly depleted in dust, though opinions strongly differ 
on the exact ext e nsion and intensity of this depletion (see 
iLi fc Greenberel Il998l for a detailed discussion on this 
topic) . To the presen t day, o ne of the most complete stud- 
ies remains that of iLi fc Gr cenbcre ( 1 99&) . taking into 
account a large set of parameters and especially realis- 
tic grain properties (porosity, size distributions, chemical 
compositions) based on observations, laboratory experi- 
ments and dust collection into space. This work claims 
that there is no more than 6 10^^ g of dust in the r < 40 AU 
region, as compared to 6 lO'^^ g in the [40,100] AU area. 
The main problem is that such SED fits are strongly model 
dependent, and in particular that the dust surface density 
distribution cannot be uniquely determined because of its 
coupling to the grain size dis tribution and to t he op tical 
properties. Furthermore, the iLi fc Greenberd l|l998l) fit- 
ting has been performed assuming that all dust is of pure 
comet-evaporation origin and that its size distribution fits 
in situ dust observations around the Halley comet. As will 
be discussed later on (section 5), this assumption probably 
cannot hold for the inner Beta-Pic disc. 

There is nevertheless one independent evidence for an 
inner dust d e pletio n deduced from direct observations: 
IPantin et alJ l|l997l) obtained resolved 12/im images of 
the r < 100 AU region, with a resolution of ^ 5 AU af- 
ter deconvolution. They concluded that there is a density 
drop of almost an order of magnitude in the innermost 
r < 10 AU area, although a puzzling density peak seems to 
be observed at 5 AU. These authors inferred a total dust 
mass of ~ 2.4 10^^ g for the r < 10 AU area. Note that 
this estimate is also strongly model dependant, though 
it doesn't make any assumption concerning the mecha- 
nism producing the dust: the authors suppose a power 
law for the size distribution with a change of power law 
in dex at a g i ven si ze (and thus 3 free parameters). The 
iPantin et al.l mass estimate is si gnificantly lower 

than the one that can be deduced from lLi fc Greenberd 



l)1998l) for the same region, i.e. ~ 2.5 10^^ g, especially 
when taking into account the fact that the upper grain 
size limit of iLi fc Gr eenl;)^r^ lll998l) 0.4 mm, is smaller 
than the 1 mm limit ofTan tiii et al.l l)l9'97l) . Extrapolating 
the Li fc Greenbera (,,1998,) estimate up to the 1 mm limit 
leads to a total dust mass of ~ 3.5 10^^ g. But as men- 
tioned before, all authors agree on one core assumption: 
there is a dust depletion in the inner /3 Pictoris system. 
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Table 1. Summary of mass estimates for the inner 10 AU region, as derived from previous works 

Authors Modeling frame Size range Mass 



Li & Greenberg(1998) full SEP fitting 



[O.l/im, 0.4 mm] 



2.5 10^^ g 



when extrapolated with a R law 



[10, 50] km 



2.510"^Mf, 



Pantin et al. (1997) 



inversion of mid-IR 
surface brightness profile 



[O.lfim, 1 mm] 



when extrapolated with a R law 



[10, 50] km 



2.4 10^1 g 



210"''Mf, 



Thebauh & Beust (2001) FEB scenario 



[10, 50] km 



15 - 50 M(f 



1.3. Parent bodies in the inner disc 

As pointed out by lArtvmowicd l)l99l1) for instance, the 
dust is not primordial and must be constantly replen- 
ished by larger bodies. This author favours the scenario of 
dust production through coUisional erosion rather than 
cometary evaporation, an assumption that we also be- 
lieve to be the most reasonable one for the inner disc 
(see discussion in section 5.3). It is then tempting to use 
this link as an indirect way to constrain the population 
of parent planetesimals, which is otherwise totally invis- 
ible to observations. The most simple way to do this is 
to suppose that a coUisional equilibrium dN oc R~^-^dR 
power law applies from the parent bodies down to the ob- 
served micron-sized grains. This is the usual assumption 
commonly made to easily derive parent bodies m asses in 
extrasolar dust discs (e.g. lAueereau et alJ l|l999^ for HR 
7496 A). This would here lead to a mass of objects in the 
10 to 50 km range compr ised between 21O~'^M0 (tak- 
ing the lPantin et al.Ul997l dust d ensity) and 2.5 10^^ 
(taking the lLi fc Greenberglll998t estimate). These values 
seem very low, especially compared to the 14 Mm mass of 
plane tesimals estimate for the whole system ijArtvmowicj. 
1997^ which was in accordance with the picture of an 
" early Solar System" . This tends to reinforce the image 
of a strong mass depletion in the inner disc. Note that 
the concurrent cometary evaporation scenario also gives 
a quantitative link between the observed du st and the 
sourc e kilometre-sized comets (e.g. Equ.2 of iLecavelieil 
1998() . but this estimate does not constrain the number of 
non-evaporating objects. 

There is nevertheless another way to get informations 
on the planetesimal population through the study of the 
so called Falling Evaporating Bodies (hereafter FEB) phe- 
nomenon. It is indeed believed that the evaporation of at 
least kilometer-sized bodies is responsible for the transient 
absorption features regularly observ ed in various spectral 
lines: Call, Mgll. Fell, etc... (e.g. iBoggess et all Il99ll 
IVidal-Madiar et alJ.ll99llBeust et al.Lll996|) . Several the- 
oretical and numerical studies have shown that these FEB 
might be bodies located at the 3:1 and/or 4:1 resonances 
with a giant planet on a slightly eccentric orbit located 
around 10 A.U. These objects are excited on high ec- 



centricity e orbits which allow them to pass sufficiently 
close to th e star, i.e. less than 0.4 AU , for silicate to evap- 
orate (see iBeust fc Morbidelll l2000t iThebault fc BeuslJ . 
l200li and references therein). iThebault fc BeustI l)2001r) 
estimated that the number density of planetesimals re- 
quired to fit the observed rate of absorption features would 
lead to a mass of ~ 15-50 objects in the 10 to 50 km 
range when assuming an equilibrium differential law in 
jj-^ ^Y^G i nner <10 AU region. This very high estimate 
is close to the Artvmo wicd ()l997f) estimate for the whole 
disc and strongly exceeds, by at least a factor lO'^, from 
the above-mentioned much lower dust-mass-extrapolated 
estimations for the inner disc. 

1.4. The need for a numerical approach 

There is thus yet no coherent picture of the inner /?-Pic 
system's structure, especially for the crucial link between 
the observed dust and unseen bigger parent bodies. The 
main reason for this is that deriving mass estimate from a 
simple power law from the micron to the kilometre might 
be strongly misleading. 

A first argument is that a very small difference in the 
power law index leads to enormous differences when ex- 
trapolating it over such a wide size range. If q is this index, 
then the mass ratio between 2 populations of sizes Ri and 
i?2 reads M1/M2 = [Ri / Ri)^''^^^ . As an example, the 
incompatibility between the 15 — 50 FEB mass esti- 
mate and the 2 10~'^-2.5 10~^ extrapolated from the 
observed dust density might be solved when changing the 
q index in the later extrapolation from -3.5 to -3.2. But 
the single power law approach raises also other problems. 
Firstly, there is no reason why the upper size limit of the 
coUisional cascade should be 10 or 50 km. Extrapolating 
a, q = —3.5 power law up to, say, 1000 km instead of 
50 km, would lead to a 4.5 times superior mass of large 
objects, thus reducing the magnitude of the inner mass 
depletion. But then, taking the same upper limit for the 
rest of disc (i.e. outside 10 AU) would increase the to- 
tal mass of the system to unrealistically high values (more 
than 1000 M0). In this case two different size distributions 
should hold for the inner and the outer systems. 
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Secondly, it is almost certain that a single -3.5 equilib- 
rium power law cannot hold over such an extremely large 
size range. For such a power law to apply, all particles in 
the system should have reached mutual coUisional equi- 
librium. This is perhaps not the case here, especially for 
the bigger bodies which, depending on their number den- 
sity, might have low collision rates. Furthermore, such a 
power law is theoretically achieve d onlv for an infinitely 
small lower size cutoff. As shown bv lCamoo Basratin et al.l 
()l994ll . any finite size cutoff will give rise to wavy size 
distribution structures which can strongly differ from the 
theoretical -3.5 slope. The reason for such size distribu- 
tion waves is simple: the smallest particles will be over- 
abundant since they have no smaller bodies to destroy 
them. This over-abundance will give rise to an under- 
abundance for all bigger bodies that might be coUisionally 
fragmented by these minimum-sized objects. This will in 
turn lead to an over-abundance of bigger bodies, and so 
on. This point is of great importance here since there is an 
obvious size cutoff for our system, i.e. the smallest grains 
that are not blown away by the star's radiation pressure 
and which are typically micron-sized. Apart from this cut- 
off effect, the smallest grains are also expected to have a 
very peculiar behaviour: even if radiation pressure is not 
able to remove them, it should nevertheless place them 
on highly eccentric orbits, thus augmenting their impact 
velocities and shattering power, but at the same time re- 
ducing their density in the inner region since they will 
spend most of their orbits very far away from the star. 
The physical link between dust and planetesimals is thus 
a complex one, that cannot be handled by simple analyt- 
ical power laws. 

We propose here to address these problems by perform- 
ing accurate numerical simulations. A statistical particle- 
in-a-box code is used to study the mutually coupled col- 
lisional evolution of a swarm of objects ranging in size 
from large planetesimals down to the smallest micron- 
sized grains. The code is similar to the ones developed 
for asteroid populations studies but stretches down to very 
small dust particles and takes into account the peculiar dy- 
namical evolution of micron-sized grains submitted to the 
star's radiation pressure. We detail in section |21 and 01 the 
numerical approach we use to derive a realistic grain size 
distribution at a distance of a few AU resulting from a col- 
lisional cascade in the radiative environment of /3 Pictoris. 
We explore the impact of the free parameters, along with 
the two extreme surface densities independently deduced 
from dust and gas observations of /? Pictoris, on the final 
size distribution after lOMyr (section^) We then discuss 
in section the implications of our approach and show 
how it helps to go towards a coherent view of the inner 
/? Pictoris disc (sectional). 

2. Numerical procedure 

We will here follow the classical particle in a box approach 
used bv models st udving the asteroid belt size distribution 
fe.g. lPetit fc Farinella>.1993ilV We consider a typical annu- 



lus of material in the inner disc, of radius 1 AU and located 
at 5 AU from the star. The system is divided into n boxes 
accounting for each particle size within the annulus. The 
size increment between two adjacent bins is 2^/^ . At each 
time step the evolution of the number dN^ of bodies of 
size Rk is given by 

n n 

dNk = ^ ni^j^kPij NiNjdt - ^ 'yi,kPi,k NiN^dt (1) 
ij=i 1=1 

where is the number of fragments injected into 

the k bin by an impact between 2 bodies in the i and j 
bins and pi_j the impact rate for a pair belonging to the 
same two bins. The last term of the equation accounts for 
the loss of k objects due to destructive collisions (7;_fc — 1 
for a catastrophic fragmentation impact and < < 1 
(satisfying the mass conservation condition) for a crater- 
ing event) . The temporal evolution of is then computed 
using a first order eulerian code with a variable time step. 
The key information needed for estimating riij^k, li,k and 
Pij is the dynamical state of the system, which can be 
parameterised by the average impact velocities {dv). 

2.1. Dynamical state of the system 

As stated in the previous section, there is yet no clear 
picture of the system we intend to study. We have in par- 
ticular no precise idea of the dynamical state of the in- 
ner disc. There is nevertheless some indirect evidence of 
the dynamical state in the outer parts, given by the ob- 
served thickness of the disc. The disc as pect ratio in the 
100 A U region is believed to be ~ 0.1 ijAueereau et all 
'2001'). Thus, a first order approximation of the average 
inclination of the observed dust particles would be half 
this value, i.e. ~0.05rd. However, this value only gives 
very partial information, and this for several reasons: 

1. It is not straightforward to extrapolate it to the inner 
regions of the disc, where the dyna mical conditions 
could be completely different. Indeed llCalas fc JewittI 
(|l995[l seem to observe a significant departure from 
constant disc opening for r < 60 AU, but such deter- 
minations should be taken with great care, since mea- 
sures of the disc's thickness become very uncertain for 
these inner regions. 

2. This value holds for the observed micron to millimetre- 
sized grains population. It is not at all certain that 
bigger parent bodies have the same inclinations 

3. The dynamical state of the system depends on the in- 
clination and eccentricity distributions. The (e) value 
cannot be directly deduced from (i), at least for the 
micron-sized population, where orbits might strongly 
depart from the equilibrium (i) = (e)/2 equipartition 
relation (points 2 and 3 will be discussed in more de- 
tails in section inj. 

We will thus take the (e) and (i) values of the parent 
bodies in the considered r < 10 AU region as free parame- 
ters (see section but we will nevertheless refer to the 
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(i) = 0.05 — (e)/2 case as our "nominal" case. Note that 
our simulations do not directly use the (e) and {i) parame- 
ter s but the average relative v elocity parameter {dv) given 
bv lLissauer fc Steward l)l993(l : 



{dv) 



1/2 



y) + {^')] {Vkep) 



(2) 



where (vkep) is the average Keplerian velocity of the bod- 
ies. Furthermore, we will assume the same (e) = (eo), 
(i) = (io) and thus (dvij) = {dvo) for all particles, with 
the important exception of the micron-sized grains which 
are significantly affected by the star's radiation pressure 
(see sectionOfor more details). We will also make the sim- 
plifying assumption that our particle in a box system is 
not dynamically evolving, so that {dv) remains constant 
throughout the run. 

2.2. Density of objects 

As described in sections 1.2 and 1.3, there are two inde- 
pendent estimates of the density of bodies in the inner 
disc: 1) a dust mass (all bodies smaller than 1 mm) of 
2.4 10^1 - 3.5 1022g derived from fits of the observed SED 
2) a mass of 15 — 50Mq of planetesimals in the 10-50 km 
range required to sustain the FEB activity. As previously 
discussed, these 2 estimates appear totally incompatible 
when assuming an equilibrium differential R~^'^ size dis- 
tribution throughout the system, since in this case the 
mass of planetesimals extrapolated from the dust estimate 
is only 2 10~^ — 2.5 10^^ Mq. In order to check how strong 
an incompatibility there really is, or if there is any incom- 
patibility at all, we will consider two extreme initial discs 
(see sections ETI and : 

— An initial R~^'^ distribution extending from 
Rmin — 2^^/^Rpr (i.c. 2 boxcs under the Rpr 
ejection size, see sectional) up to Rmax — 50 km which 
is compatible with observed dust mass estimates. We 
chose to take an initial dust mass of 2 10^^ g in the 
whole inner disc (i.e. ~ 2 lO^^g in the considered 1 
AU annulus around 5 AU), a n int ermediate value 
betwe en the IPantin et al.l (^223) iLi fc Greenberd 
l|l998 ^ estimates. This leads to a total initial mass of 
1.8 10^^ g for the whole system. This initial density 
distribution will be taken for our "nominal" case 

— An initial R^^'^ distribution extending from 
Rmin = 2^'^/^Rpr up to R^ax = 50 km and 
compatible with an average FEB estimate of ~ 25M0 
of 10-50 km-sized objects in the 1-10 AU region (i.e. 
~ 2.5 M0 in the considered annulus) , which leads to 
a total initial mass of 3.5 10^* g 



Regardless of the initial density distribution, the number 
of size boxes considered is always the same for a given Rpr^ 
ranging from Rmin = 2~^/^Rpr to Rmax = 50km, with 
two adjacent boxes separated by a factor 2^/^ in size, thus 



leading to a total number of 103 boxes for the "nominal" 
case where Rpr = 5 10~^cm (see sectional 



2.3. Threshold specific energy 

The core of such a code is the prescription giving rii.j^fc for 
a given {dv). Basically, impacts can be divided into two 
categories: catastrophic fragmentation and cratering, de- 
pending on the value of the impacting energy as compared 
to Q*, the threshold specific energy of the bodies, which 
represents their resistance to shattering and is deduced 
from laboratory experiences and analytical considerations. 

is by definition the value of the specific energy Q (the 
ratio of the projectile kinetic energy to the target mass) 
when the mass of the largest remaining fragment is 
equal to 0.5 Mi. 

The problem is that estimations of Q.^, do strongly 
differ from one author to another (see figure 8 of 
iBenz fc AsDhauglll999l for an overview). Basically, all au- 
thors agree on one core assumption, i.e. the response of 
solid bodies to impacts is divided in two distinct regimes: 
the strength regime for small bodies, where the object's 
resistance decreases with size, and the gravity regime for 
larger objects where resistance incre ases with size because 
of th e object's self-gravity (e.g. iHousen fc Holsapnlel 
Nevertheless, the slopes and turn over size from 
one regime to another are still a great subject of debate. 
We will here consider separately two different prescrip- 
tions (see section E3)l : 

— the global strength + gravity regime law given in Equ.6 
of Benz fc Asphaue ( l^ftOi^ fnom inal case) 

- the IHousen fc Holsapplel Il990tl law for t he strength 
regime completed by the lHolsappl3 ( ^3) l^^w for the 
gravity regime. 



Our code calculates for every target-impactor couple 
(i, j) the corresponding value of . Let us term Fi f{ij) 

the ratio Mif(i)/Mi. From the value of Q* , Ffia can be 
inferred through the empirical relation ijFuiiwara et alJ . 
I1977I) : 



F, 



0.5 



Er 



1.24 



(3) 



where Erei is the relative kinetic energy of the system 
given by Erei = MiMjdv^/2{Mi + Mj). Note that this 
relation is valid only for head-on impacts and has to be 
corrected by taking into account its va lue averaged over 
all impacts angles. We will here follow IPetit fc Farinellal 
(1993) and take the average value: 



F 



= 3F. 



2/3 



2F 



(4) 



Thus, Equ|31 has to be corrected by a numerical factor 
Xcr = 4-1/1-24 ^ 0.327, since T]J^ = 1/2 for 
1/8. 
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2.4. Fragmentation 

Catastrophic fragmentation occurs by definition when 
Fifi^ij) is less than 0.5. If we suppose that the produced 
fragment size distribution follows a single-exponent power 
law dN = CR'^dR, then there is a unique set of values for 
q and C derived from the value of Fif^i^j) and the mass 
conservation condition. As pointed out in several previous 
studies, this single power law specification is the easiest 
to handle in models but it is a strong oversimplification. 
It gives rise to several problems, in particular the possi- 
bility to get so-called " supercatastrophic" impacts where 
q < —4, for which there is a divergence of the total mass 
when taking infini tely small lower cutoff. As noted by 
iTanga et al.l l(l999f) "...values beyond -3 for the exponent 
of the cumulative size distribution cannot hold down to 
very small sizes, because this would lead to unreasonably 
large reconstructed masses. For this reason it is clear that, 
at some value of the size, the distributions are expected to 
have a definite change of slope" . Note that this change of 
slope between the small and large fragment s domain is also 
supp orted by experimental experiments (jPavis fc Rvari 
This problem is particularly crucial for the present 
study since our size cutoff is extremely small (see below). 

As a consequence, we will here adopt 2 different power 
laws of index qi and q2, each holding for a different mass 
range and always taken such as the small mass index q2 is 
smaller than qi . The main problem is to determine where 
the change of slope occurs and what the difference in slope 
is. We shall remain careful and keep the slope changing 



2.5. Cra taring 

For the cratering c ase (Fif > 0.5), we will t ake the simpli- 
fied prescription of lPetit fc Farinellal l(l99,'^ . where a fixed 
power law index qc = —3.4 is considered. The total mass 
of craterized mass is given by 



Mr, 



aErei for 

QXcrOt 
100(5*Q! — Xc 



rel 



< 



-Erel + 



100a ' 

Mi Xcr — 10Q*a 



10 Xrr - loos'. 



(9) 

(10) 



for : Erel > 



M, 
100a 



where a is the crater excavation coefficient which depends 
on the material properties. We will explore values of a 
fsection f4.6|) ranging from 10^^ to 4.10~* s'^.cm''^, the ex- 
treme values corr esponding to " hard" and " soft" material 
respectively (see IPetit fc Farinellal Il993l and references 
therein), and take 10~® s^.cm"^ as our "nominal" value. 
The mass of the largest fragment produced by the impact 
is then equal to Fi^^ Mcr, where Fi^^ = 1 + ^{qc + 1). 

2.6. Fragment reaccumulation 

The fraction of fragmented material reaccumulated onto 
the parent bodies is the result of the competing ejectas' 
kinetic energy and the parent bodies' gravitational po- 
tential. We will make the simplified assumption that all 
fragments produced after an (i, j) impact have the same 



io. vv^ oiicuii i^iiiciiii v.cj,iv.iui ciii^i ivv^^jj v^L^ oi^^j^^ i^iit^iigiiig, tragments produced alter an imp act na 
size as well as the ratio qu/q2^ as free parameters velocity distribution llStern fc Colwellll997i) : 



that will be explored in the runs (see section • Note 
that once Rs and qii/q2i are given, the values qu and q2i 
for the fragments produced on a target i by an impactor 
j are uniquely determined through the set of relations: 



^ fke-^rel(^i j^ 

' M 



1/2 



(11) 



Ml, = 



(1 - bu) y 'f^') 



M 



(l-bu) 
s(^) 



M, 



Cii - ibuR^fll^ 



Co 



3buR 



.3bii 



n3(bi.-fc2.) 



C2i 



(3 - 362O 



R 



-362 



M, - Mu 
Ms{i) 



(5) 



(6) 



(7) 



(8) 



where fke is the fraction of kinetic energy that is not 
dissipated after an impact. We will make here the clas- 
sical assumption that fke = 0.1 for high velocity impacts 
l|Fuiiwara et alJ . Fl989|) . The mass fraction of fragment ma- 
terial that escape s the t arget -|-impactor system is given by 
llStern fc Colwelll. Il997l) : 



fe. 



0.5 



Vfr 



-1.5 



(12) 



where bu = -|(gii + 1) and 62; = -3(921 + 1), ^^(i) is 
the mass of an object of size Rs(i), Mu is the total mass 
of fragments produced in the domain where the dN — 
CiiR'^^^dR law applies, i.e. between the size of the slope 
transition Rs(i) and the size of the largest fragment Rif^i^, 
and Cii and C2i are the coefficients for the two dN = 
CuR'^^'dR and dN = C2iR''^'dR power laws. This set of 
equation is solved numerically for each (i, j) couple. 



where Vesc is the escape velocity of the colliding bodies 
system. 

3. The specific behaviour of the micron-sized 
population 

The main challenge of this simulation is that we would like 
to study the coUisional correlation between objects rang- 
ing from the micron-sized to the kilometre-sized domain, 
i.e. separated by a 8 orders of magnitude in size. Our lower 
cutoff is indeed the " real" physical cutoff Rpr imposed by 
the effect of the star's radiation pressure. Radiation pres- 
sure also strongly affect particles bigger than i?pr, but 
still in the same size range, by placing them on highly 
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eccentric orbits. These eccentricities depend on the ratio 
(3 between the radiation pressure force Fpr and the grav- 
itational force Fgrav ■ For a particle produced by a parent 
body on a (cq, Cq) orbit at a distance rg from the star, one 
gets : 



1-/3 



1 - 2ao/3/ro 



ao 



(l-2ao/3/ro)(l-e§) 
(1-/3)2 



1/2 



(13) 



(14) 



where Ofc and Ck are the produced grain's semi-major axis 
and eccentricity and ao and eg the semi-major axis and ec- 
centricity of the parent body. From Eauation ll4l and with 
the assumption that the planetesimals releasing dust par- 
ticles by collisions are mostly on circular orbits (eo — 0, 
Oo — ro), then grains with /3 > 0.5 are ejected from 
the system on hyperbolic orbits. The cutoff size Rpr is 
by definition the size for which /? = 0.5 and grains with 
sizes R larger than Rp^ are related to /3 through the re- 
lation (3 = 0.5{Rpr/ R). The blow-out size Rp^ depends 
on the stellar spectra and on grains optical properties. 
For /3 Pictoris we use a low-resolution A5V spectra de- 
rived from Kurucz stellar models and we adopt the chem- 
ical grain composition proposed by[ A fc Greenberd lIlQQS ) 
(we refer to the latter paper and to Aiigereau et all Il99f , 
for a full discussion on chemical and optical properties 
of the grains assumed here) . Bare compact silicate ( "Si" ) 
grains in the surroundings of /3 Pictoris and smaller than 
-Rpr, compact — 3.5 /im havc /3 > 0.5. The same grains but 
coated by an organic refractory ( "or") mantle in a Si:or 
volume ratio of 1:2 as proposed bv lLi fc Greenber3 l)l998|) 
are ejected from the system on unbound orbits if they are 
smaller than i?pr, compact — 2.5 /im. Actually the main un- 
certainty on Rpr relies on the grain porosity P. The poros- 
ity affects the optical properties of the grains and conse- 
quently Frad- But actually the /3 ratio more dramatically 
depends on P through the grain density in Fgrav especially 
for large porosities. The density of porous grains is re- 
lated to the density of the same but compact grain by the 
simple relation : Pporous = (1 - -P)/Ocompact which implies 



for large porositie s : R 



{l-P)-^R, 



^ ^ I -n-r norniis — V J- i_J ■-pi-, compact ■ 

From SED fitting, iLi &: Greenberd l|l998|) constrained P 
in the narrow range [0.95, 0.975]. But these values are 
obtained when assuming that all grains are of cometary 
origin, an assumption that we believe might not hold for 
inner /3 Pictoris disc (see the more complete discussion in 
section 5.3.) In the present paper we keep P has a free 
parameter with Rpr — 5 /im (i.e. P ~ 0.5) taken as our 
reference nominal case (see section 

The radiation pressure induced eccentricity expressed 
by Equation 1141 is significant, say > 0.1, for all particles 
comprised between Rpr and ^ bRpr- Thus all objects in 
this size range will have orbital characteristics that depart 
from the general average values defined in section 2.1. This 
will significantly affect the collision rates and physical out- 
comes for impacts involving these small grains. For these 



Table 2. numerical estimate of {fnk)) in the 5 AU region, 
for a swarm of grains produced in the whole 1-10 AU re- 
gion, as a function of [3k (see text for details). Note that 
for low /3, Ck might become lower than the average eccen- 
tricity for the parent bodies in the disc (see section 2.1). 
In such a case, the radiation pressure effect is neglected 
for all Ck < (eo) and the values of (fi^k)) rescaled so that 
iliik)) = 1 for Ck = (eo). 







{fi{k)) 


0.49 


0.96 


0.038 


0.39 


0.63 


0.344 


0.31 


0.45 


0.521 


0.24 


0.32 


0.665 


0.19 


0.24 


0.753 


0.15 


0.18 


0.812 


0.12 


0.14 


0.849 


0.10 


0.11 


0.887 


0.08 


0.09 


0.930 



impacts, Equ|2] is no longer valid in its simple form and 
(dvij) will be numerically estimated. To do this, we use 
a simplified ver s ion of a deterministic coUisional model 
l)Thebault et al.L l2002l and references therein) to derive 
average (dv) between a population of targets having the 
nominal orbital characteristics as defined in section 2.1 
and a population of impactors with a given /3fc (i.e. ak and 
efc), all randomly distributed within the 1-10 AU region. 

Another major consequence of these radiation-induced 
high efc is that small grains will spend a significant fraction 
of their orbits outside the inner disc. Thus, at a given mo- 
ment, only a fraction fi(^k) Nk of these bodies will actually 
be present in the considered system. These {fi{k)) are nu- 
merically estimated with a simple code randomly spread- 
ing 10000 test particles of a given Pk uniformly produced 
in the 1-10 AU region (table 0). 



3.1. Timescale 

It is important to note that these {fi[k)) values are not 
reached instantaneously: small grains produced after an 
impact need time to reach the remote aphelion of their 
high a and high e orbits. If dtej{k) is the typical time 
needed for a small grain produced in the inner disc to reach 
r ~ lOAU when placed on a high Uk and ek orbit, then 
during a time step dti , the fraction of produced k grains 
that leaves the system will be approximated through the 
simplified relation: 



-(l-/.w).(l-exp -^.c-)) 



(15) 



Bodies that did not leave the system during dti are then 
added to a "bodies on their way to leave" subdivision of 
the k population, that will in turn decrease by a (1 — 

/i(fc))-(l ~ exp '''^jC:) ) fraction at the next time step dt2, 
etc... 
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3.2. Collisional destruction outside ttie inner disc 

Another possible effect affecting tlie smallest high /3 par- 
ticles is that a fraction of them might be destroyed by col- 
lisions outside the considered inner disc, since they spend 
an important fraction of their orbit close to their apoas- 
tron which might lie beyond 10 AU. These collisions would 
prevent them from re-entering the inner system. Such col- 
lisions are by definition not modeled by the collisional evo- 
lution EquQ] and this could lead to an overestimation of 
the density of these bodies. 

Taking into account this apoastron-coUisions removal 
effect requires to make physical assumption on the 
external parts of the disc and would add several badly 
constrained free parameters to our already large set of 
variables, in particular the rates and timescales for these 
collisional destructions as a function of /3. We nevertheless 
tried to investigate the possible importance of this effect 
by performing test runs where we artificially introduced a 
new parameter fci{k) i standing for the fraction of k bodies 
destroyed by collisions in the r > 10 AU region and the 
corresponding destruction time scale dtci(k)- The induced 
removal of k bodies is then treated the same way as in 
EquEl fci(k) might be taken equal to the fraction of j3k 
grains produced within the inner 1-10 AU disc which have 
their periastron outside 10 AU. The dependency of dtci(k) 
with f3k is more difficult to establish, and several values 
will be explored. 

As will be shown in section 4.7, the obtained results 
do not significantly depart from our "nominal" case. As a 
consequence, and for sake of clarity, we chose to neglect, 
in a first approximation, this apoastron-coUision effect. 

3.3. Particles with /3 > 0.5 

Even if these particles' ultimate fate is to leave the system, 
their ejection takes also a certain amount of time and nu- 
merous very small grains, in this transition phase towards 
ejection, might be present in the system and thus coUi- 
sionaly interacting with other objects. As a consequence, 
our runs will be performed with 2 bins below the limiting 
(3 = 0.5 size. The fraction of /3 > 0.5 bodies that do not 
leave the system after an impact is computed the same 
way as in Eau ll5l by numerically estimating dt^j^k) and 
setting 0. 

4. Results 

We present here the results obtained for several runs ex- 
ploring all important parameters the system's collisional 
evolution depends on. As previously mentioned, we define 
as our "nominal" case the one defined in section 2.1. Initial 
conditions for this reference case are summarized in Tab|31 
For sake of clarity, all other parameters are separately 
explored in individual runs, even though some parame- 
ters should in principle not be independently explored, 
like in particular the value of Rpr (i.e. the grains' poros- 
ity) and the fragmentation and/or cratering prescriptions. 



Table 3. Initial parameters for the nominal reference run 
(see text for details) 



Minimun size bin 
Maximum size bin 
Rpr {13 = 0.5) 

(eo) 
(io) 

Qt law 
&1/62 

Rs{i) 

excavation coefficient a 
Initial density 



3.1510-*cm 

5.4 let'em 

510~''cm 

0.1 

0.05 

Benz fc Aspliaug (1999) 
1.5 

^{,)/2.10'5 
10"**s^cm~^ 

dN oc R~^-^dR distribution 

with Mdust = 2 10^^ g in the annulus 

{Mdust ~ 2 10^^ g in the whole inner disc) 



2 bOO years 
25 000 years 
500 000 years 

3 000 000 years 
10 000 000 years 




size (cm) 

Fig. 1. Size distribution for the low-mass nominal system 
(see tablflSI at 5 different epochs. Note that the y-axis dis- 
plays the mass contained in one size bin, which is a correct 
way of displaying the mass distribution since all size bins 
are equally spaced in a logarithmic scale. This plot is more 
"visual" than a more classical dN{R) one, since it can be 
directly interpreted in terms of mass contribution (and 
mass loss or mass increase) of each size range to the total 
mass. 



All runs are carried out until t final = 10^ years, i.e. ap- 

Siroximatcly one third of the minim um age of the system 
Barrado v Navascucs et al.Lll999|] . 

4.1. Nominal case 

Fig- n shows clearly how the system quickly departs from 
the initial R~^'^ distribution. A wavy structure rapidly 
appears be cause of the minimum size c ut-off, in accor- 
dance with ICampo Baeatin et all l)l994|) . This structure 
is building up progressively, starting from the lowest sizes 
and expanding towards the bigger objects bins. A quasi 
steady-state is reached after ~ 10® years and no signif- 
icant further evolution of the system is observed in the 
next 9 10® years, except for a slow decrease of the system's 
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Fig. 2. Temporal evolution of the system's total mass for 
different cases 



total density. As could be logically expected, the wavy 
structure is the most significant in the small size domain. 
There is in particular a strong mass depletion, of a factor 
simeqAO, in the 10~^-1 cm range, with the lowest density 
point around R ~ 0.1cm. This depletion has 2 distinct 
causes: 1) The overabundance of very small particles due 
to the size cut-off. Note however that this overabundance, 
though still present, is significantly damped or even erased 
for the smallest particles (close to Rpr), because these bod- 
ies spend a significant fraction of their orbits outside the 
inner disc (see the fi(k) parameter in TabEJ- 2) The high 
(dv) values for impacts involving particles close to Rpr, 
which are on highly eccentric orbits. 

Another important result is that the total mass loss 
of the system over 10^ years remains relatively limited, 
i.e. less than 12 % (cf. FiglSJ- Furthermore, the ratio 
Mdust/Mpianetesimais, after large initial variations, pro- 
gressively converges towards a value which is ~ 1/3 of 
the dN = CR~^'^dR power law value. This is mostly due 
to a decrease of Mdust, with Mpianetesimais being almost 
constant. 



Fig. 3. Temporal evolution of the ratio 

Mdust /Mpianetesimais for different cases, where Mdust 
is the total mass of all objects smaller than 1 mm and 
Mpianetesimais IS the mass of all objccts bigger than 10 km. 
The horizontal line gives the initial value corresponding 
to an academic dN = CR^^'^dR size distribution. 
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1 1 1 1 1 1 
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.1=4 000 yeors 
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. . . t = 300 000 years 


\- / / 


... , t = 3 000 000 yeors 




1 , 1 


t = 10 000 000 yeors 



10 10 ^ lO'-' 10"^ 10'* 10^ 10^ 

size (cm) 

Fig. 4. Size distribution at 5 different epochs for the 
massive-disc case, where the total mass of the system is 
chosen in order to match the planetesimal mass estimates 
deduced from FEB mechanism analysis (see section 2.2). 



4.2. Massive disc 

The massive-disc case turns out to be significantly dif- 
ferent (FigQ. Although a quasi steady-state is here also 
rapidly reached and has a profile similar to that of the 
nominal case, this steady-state is obtained for a much 
higher density. This leads to much faster mass loss than in 
the previous case, exceeding one order of magnitude at the 
end of the simulation (FigO, since mass loss in a given 
collisional system increases with the square of the system's 
density. We discuss the implications of these results on the 
FEB phenomena in section 1531 



4.3. Role of tlie disc's excitation 

Apart from the nominal case with (i) = l/2(e) — 0.05, 
two different disc excitations have been tested: one low 
excitation case at (i) = 0.0125 = l/2(e) and one high 
excitation case at (i) — 0.1 = l/2(e), all other parameters 
being equal (FiglHJl. 

As could logically be expected, the total mass loss is 
much higher in the high excitation case, ~ 18%, than in 
the low excitation case, ~ 4%. Nevertheless, the steady- 
state regime profile is significantly different for both runs. 
The density well in the 0.01-1 cm range is in particular 
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size (cm) 



Fig. 6. Final distribution (at t — lO^years) for different 
levels of the disc's dynamical excitation 



much deeper for the dynamically cold disc. This is a fully 
logical result when considering the fact that the radia- 
tion pressure induced high e of the smallest grains only 
weakly depends on the dynamical state of the parent bod- 
ies ('Eau ll4|l . As a consequence, the contrast between the 
excitation, and thus the shattering power, of the smallest 
grains and that of the rest of the particles is very high. 
The destruction rate of bodies in the 0.01-1 cm range by 
small grains is thus at the same level than in the nominal 
case, whereas the production rate of 0.01-1 cm grains by 
collisions between bigger objects is much lower, hence the 
deeper density well. Conversely, for the high excitation 
case, the contrast between the small grains' and bigger 
objects' destructive powers is significantly damped, hence 
a shallower density drop in the 0.01-lcm range. 



4.4. Porosity, value of Rpr 

As previously discussed, our nominal case corresponds to 
compact low porosity grains and Rpr = 5 fim. We inves- 
tigated different porosities, and thus different Rpr values, 
all other parameters being equal (with always 2 size boxes 
below Rpr). As can be clearly seen in Figd changing the 
value of Rpr results mainly in shifting the wavy struc- 
ture without affecting its overall profile. Although it is not 
strictly speaking a homothetic shift, mainly because of the 
complexity of the prescription, differences are minor 
ones, and the density drop is always located at roughly 

100 Rpr. 

4.5. prescription 

Taking the Housen & Holsapple (1990) and Holsapple 
(1994) prescription for leads to a final size distribution 
which is remarkably close to the nominal Benz & Asphaug 
(1999) case (FiglSjl. The main difference is a more defined 
density drop for objects bigger than 10^ cm, which is di- 
rectly due to the fact that values for bodies in this size 
range are significantly lower than in the Benz & Asphaug 
(1999) model. This deficit in large bodies is the reason for 
the more significant total mass loss in the system (Figl^J, 
since these bodies contain most of the system's mass. Note 
however that this total mass loss does not reach very large 
values, remaining limited to less than 30%, and that the 
global dust to planetesimals mass ratio is also relatively 
close to the nominal case value (FigOJ. 

4.6. Fragmentation and Cratering prescriptions 

As previously shown, the main free parameters for our 
fragmentation prescription are the ratio qi / q2 and the size 
Rs of the slope transition. These parameters were both 
explored in independent runs whose results are presented 
in Figjnj. As can be clearly seen, the values of qi/q2 and 
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lO'sl i , \ , \ , \ , \ , \ , 

10-^ 10-^ lO'' 10^ 10^ 10^ 10^ 

size (err) 

Fig. 8. Final mass distribution (at t = 10^ years) for dif- 
ferent Q* prescriptions 

Table 4. This table sums up, for all objects within 3 
different size ranges, the respective amount of mass that 
is removed by all cratering and all fragmenting impacts. 
These values are obtained in the steady state regime for 
the nominal case. 



size range (cm) 


R<5. lO-" 


0.01<R<1 


10='<R 


total 


fraction of 










mass removed: 










by fragmentation 


0.89 


0.13 


0.76 


0.67 


by cratering 


0.11 


0.87 


0.24 


0.33 



Rs only moderately affect the final size distribution within 
the system. 

As appears in Fig (|10|l . the cratering prescription, in 
particular the value of the excavation coefficient a, has a 
more significant effect on the physical evolution of the sys- 
tem. Taking a very hard material prescription {a = 10^^) 
leads indeed to a final size distribution which is very close 
to a dN — CR^^'^dR power law. The density drop in 
the sub-centimeter size-range is in particular significantly 
reduced, with only a factor 8 drop in a narrow region 
around 10^^ cm. Conversely, the very soft material run 
(a = 4.10"^) leads to a deeper density drop and a more 
pronounced wavy structure throughout the size distribu- 
tion. This dependency of the size distribution profile on 
the cratering prescription is easily understandable when 
realizing that, in the 0.01 to 1cm domain, cratering is a 
much more efficient process than fragmentation in terms 
of mass removal (Tab0J; mainly because of the cratering 
events due to the high e grains in the Rpr to ~ 10 Rpr 
range. Note however, that for the system as a whole, it 
is fragmentation which is clearly the dominating mass re- 
moving process (Tab^. 





1 1 1 1 1 1 1 1 

nominal case bi/b2=1.5 and Rs=R(i)/2. i 0^ 






b,/b2=1.15 


- 




b,/b2 = 2 

. . . R3=R(i)/2.103 








- 










1,1,1,1, 





10-4 10-2 -10° 10^ ^o'^ 10^ 10^ 

size (cm) 

Fig. 9. Final mass distribution (at t = 10^ years) for dif- 
ferent values of the free parameters of our bimodal power 
law: i.e. the ratio of their slopes qi/q2 and the position of 
the slope changing size Rs with respect to the size Ri of 
the impacted body. 



10=" 




lO-l 10-2 10° IQZ lO* 10° 10= 

size (cm) 

Fig. 10. Final mass distribution (at t = 10^ years) for 
different values of the excavation coefficient a 

4. 7. Collisions outside the inner disc 

As discussed in section 3.2, we chose to perform additional 
test runs checking the possible influence of small particles 
removal by collisions outside the inner disc. This effect 
is arbitrarily parameterised by the two parameters dtci{k) 
and fci{k) (see section 3.2). 

We present here results obtained for the most extreme 
case, where dt^n^f^-^ was unrealistically supposed to be equal 
to one orbital period of a j3k particle. As appears clearly 
from Fig llll differences with the nominal case remain 
marginal. As expected, the main difference is found for 
bins just below the (3 — 0.5 cut-off, with a factor 4 num- 
ber density difference for the first bin corresponding to 
bound orbits (/3fc = 0.49). Nevertheless, this difference al- 
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xtreme opoastron-collisions effect 




10^ 
size (c 



Fig. 11. Final mass distribution (at t = 10^ years) for 
an academic case with a very efBcient removal of small 
particles by hypothetic collisions in the r > 10 AU region 
(see text for details). 

ready drops to 25% for particles of size 2Rpr {Pk — 0.24). 
As a consequence, the sharp density drop induced by the 
apoastron collisions effect remains confined to a narrow 
size range of particles. Furthermore, this density drop does 
not have significant consequences on the rest of the size 
distribution and doesn't affect the global profile of the 
wavy size distribution. This is because it affects parti- 
cles which are already strongly depleted because of their 
high a and e (low /i(fc)) values. Besides, this removing ef- 
fect's dependency on f3 is relatively similar to that of the 
one induced by low f^k)) values; it will thus only tend to 
reinforce an effect already taken into account. Thus, fur- 
ther depleting these populations does not lead to drastic 
changes. 

5. Discussion 

5.1. The massive disc case. Problems with the FEB 
scenario 

One of the most obvious and easily understandable re- 
sults is the strong mass loss for the massive disc case. 
As previously stated, the system loses more than 90% 
of its total mass over 10^ years. The problem gets even 
worse when trying to extrapolate this mass loss to the 
past. This might be done when noticing that, apart from 
the initial transition phase, the system's mass loss in the 
steady state regime might be fitted by a M(() = (a + b.t)^^ 
law (which is a logical result since the mass loss is pro- 
portional to the square of the system's total mass), with 
a = 5.210^2^5^1 and b = 3.210^35^-1^-1^ j^- jg g^gy 

to see that extrapolating this law to the past leads to 
masses that become rapidly unrealistically high. Thus, the 
planetesimal density required to sustain the FEB phe- 
nomenon corresponds to a very rapidly evolving system 
and cannot be maintained over a long period of time. 



One could argue that considering a dynamically colder 
system would significantly reduce the system's mass loss. 
But this would not solve the problem, since the strength 
of the FEB producing mean-motion resonances directly 
depends on the system's excitation, so that reducing the 
disc's excitation would require an even higher number den- 
sity of planetesimals in ord er to get the observed FEB rate 
l|Thebault k Beustl l200l|) . 

The problem with the sofar accepted FEB scenario is 
then the following: from combined observations and mod- 
elling, the massive disc required to sustain the observed 
activity should erode significantly within less than lO^yrs, 
giving a natural end to the FEB phenomenon. If this was 
to be the case, then we would be presently witnessing a 
very transient phenomenon. This does not appear satis- 
factory from a statistical point of view. Should this mean 
that the FEB scenario should be rejected as a whole ? 
We believe that it is too early to state anything definitely. 
There are several reasons for that: 

We must first recall that the estimate for the neces- 
sary disc population for sustaining the FEB activity is 
derived through a chain calculation which depends on sev- 
eral poorly constrained parameters (s e e the extended dis- 
cussio n in iThebault fc Beustl l200l|) . iThebault fc Beusd 
(I2OOII) (Eqs. 7 and 8) showed that the most crucial pa- 
rameter is here i?FEB, i-e. the minimum size of bodies 
able to become observable FEBs, since the deduced total 
mass scales roughly as i?pEB ™ simplified case where a 
power law of index q applies for the size distribution above 
-Rfeb, so any change to R-p v.-n may induce drastic cha. nges 
to the estimated disc mass. IThebault fc Beustl l)200lll as- 
sumed i?FEB = 15 km, but this value is poorly known and 
could easily vary by one order of magnitude. i?FEB exists 
because bodies smaller than i?FEB are assumed to evapo- 
rate too quickly and consequently make too few periastron 
passages in the refractory evaporation zone 0.4 AU) to 
significantly contribute to the observable spectral activity. 
The value of -Rfeb is thus related to the evaporation rate 
of the FEBs themselves. Simulations of the dynamics of 
the material produced by FEB evaporation (Beust et al. 
1996) led to derive production rates of a few lO^kgs"^ as 
necessary to yield observable spectral components. We be- 
lieve that this part of the scenario needs to be revised. The 
main reason for that is that in Beust et al. (1996) simu- 
lations, the material escaped from the FEBs was assumed 
for simplicity to consist of the metallic ions we study and 
volatile material. The metallic ions undergo a strong radi- 
ation pressure from the star while this is not the case for 
the volatiles. Hence the volatiles retain the metallic ions 
around the FEB coma for a while, leading to observable 
components. The production rate was then derived from 
the necessary amount of volatiles to retain the ions, and 
from assumptions about the chemical composition of the 
body. All t his is obviously th e weakest part of the scenario. 

Besides. lKarmann etaP (JOOl) showed that if the FEB 
progenitors are supposed to originate from 4-5 AU from 
the star, they should no longer contain ices today (i.e. 
volatile material), apart from an eventual residual core. 
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More recently, iKarmann et al.l l|2003r) made an indepen- 
dent theoretical study of the evaporation behaviour of 
such objects when they gradually approach the star on re- 
peated periastron passages. The evaporation rates derived 
are thus independent from any observation. Basically, this 
work shows that 10 km sized bodies fully evaporate with 
repeated periastron passages, and that evaporation rates 
of a few 10''kgs~^ are actually reached, but this occurs 
only when the periastron is less than ~ 0.2 AU, i.e. well 
inside the dust evaporation zone and shortly before the 
final evaporation of the body. Before that, any FEB en- 
tering the dust evaporation zone 0.4 AU) but for which 
the periastron has not yet reached 0.2 AU actually evapo- 
rates, but at a weaker rate. If it is small, it thus survives 
more periastron passages than in previous estimates, and 
may contribute to the observational statistics. 

However, whether bodies with no or very few volatiles 
may generate observable components is questionable, as 
volatiles have a crucial role in the dynamics of the metallic 
ions. Within the refractory material, some species suffer- 
ing low radiation pressure, and that are probably abun- 
dant (Carbon, Silicon,. . . ) may play the retaining role of 
volatiles. Obviously this question must be reinvestigated 
with more realistic simulations, but a probably outcome 
will be that -Rfeb could end up to be at least one order of 
magnitude less than previously estimated. In this context, 
our chain calculation would lead to a much lower disc mass 
necessary for sustaining the FEB activity. 

In this context, it is impossible to rule out the FEB 
scenario on this basis alone, but this remain a problem- 
atic possibility. All we can presently state i s that the disc 
populations inferred bv lThebault &: Beusd l)200l|) are un- 
realistic and that the FEB scenario should at least be 
reinvestigated much more carefully. 

5.2. Departure from the R^'^'^ profile 

Putting aside the peculiar massive disc problem, the most 
striking result, present for almost all tested simulations, is 
a final size-distribution that significantly departs from a 
R~^'^ power law, especially in the small size domain. The 
only exception to this behaviour is a run with a very hard 
material parameter for the cratering prescription, which 
means that alternative size-distribution profiles cannot be 
completely ruled out, although they seem to represent a 
marginal possibility. Of course, due to the complexity of 
the studied problem, all free parameters could not be ex- 
haustively explored. Besides, there are some parameters 
that are strongly coupled, i.e. fragmentation and crater- 
ing prescriptions should in principle not be independently 
explored. Nevertheless, there seems to be a global ten- 
dency towards a common feature which consists of a lack 
of objects close to the limiting ejection size Rpr, a density 
peak at ~ 2 Rpr and a sharp density drop compared to 
the i?^"^-^ law, of one or two orders of magnitude in mass, 
at ~ 100 Rpr. It is also important to note that a very 
badly constrained parameter such as the disc's dynamical 



excitation does not seem to have a crucial influence on the 
profile, thus reinforcing the genericity of this result. 

These departures from the dN = CR~^'^dR power 
law do not lead to radical changes in the global dust vs. 
planetesimal mass ratio in the system, which only de- 
creases by a factor ~3-4. If 2 10^^ g is a typical value 
for the amount of dust (i.e. i? < 1 mm) in the inner 10 
AU region (see section I), then we estimate from our re- 
sults that the corresponding mass of 1 km<R< 50 km ob- 
jects should be ~ 3.5-7 10~^ M^, which remains a value 
comparable to the one roughly derived from a R~^-^ 
power law (see section I). Even stretching this value up 
to the 1 km<R< 500 km range does not lead to more than 

0. 15 Af® of "large" objects. Our calculations thus quanti- 
tatively confirm what had been previously inferred from 
questionable assumptions (an R~^'^ power law): there is 
a lack of objects, that holds even for large planetesimals, 
in the inner disc. 

This is an additional problem for the FEB scenario, 
since this value is far from being enough in order to ac- 
count for the sharp incompatibility between the amount 
of observed dust and the required amount of FEB induc- 
ing planetesimals. In any case, it appears clearly that the 
FEB model as it is currently accepted cannot be compat- 
ible with a " reasonable" estimate of the dust production 
rate in the inner disc. 

5.3. Collisional erosion vs. cometary evaporation 

Let us recall that th e precise SED fit performed by 
iLi Sz Greenberd l)l998f) was obtained assuming that the 
dust is of pure cometary origin and is not affected by col- 
lision processes. On the contrary, in our simulations we im- 
plicitly made the assumption that the inner /3 Pictoris dust 
disc is made of collisional debris. We do believe that our 
results retrospectively justify this assumption, although 
without ruling out the possible presence of evaporating 
bodies, and this for several reasons. 

1. For what is presently known about the inner disc, the 
collision production hypothesis seems to be quanti- 
tatively more generic t han the concurrent cometary 
evaporation. In his Ecu . 2 iLecavelier (1998,) proposed a 
simple expression in order to estimate the number Nco 
of currently evaporating comets in a dust disc, where 
Nco directly depends on Mdust, the total dust mass, 
and the typical lifetime for a dust particle before 
destruction. Takin g the s ame Halley-at-l-AU evapo- 
ration rates as iLecavelied (|l998) and bodies of 20 km 
in radius leads to Nco — 1-510* * (tdust/lC yrs)^^. 
Considering that collisional lifetimes of dust grains 
are of the order of 10'^ years in the inner Beta-Pic 
disc fFig ll2|l . we get ~ 2.6 10"'^ of evaporating 
comets in the r < 10 AU region. This value repre- 
sents ~ 10% of the estimated total mass of kilometre- 
sized objects (cf. Tab0, which should mean that 
one out of ten planetesimal objects is an evaporating 
comet. Although such a possibility cannot be com- 
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Fig. 12. Typical lifetime, as a function of its size, of a dust 
particle in the inner disc before destruction by collision 
(steady state regime of the nominal case) 

pletely rul ed out, it seems nev ertheless rather un- 
likely since iKarmann et al has showed that all 
kilometre-sized objects originating from the 5 AU re- 
gion should no longer contain volatile material today. 
One could argue that the previous calculation depends 
on poorly constrained parameters and that the re- 
quested number of evaporating bodies could be lower. 
But even in this case there is one crucial problem to 
solve : what is the mechanism constantly refilling the 
inner disc with so many fresh comets? The main diffi- 
culty is that this refilling has to be very effective and 
rapid, since the average time needed for a 10 km body 
at 5 AU to lose all its volatile material is less than 
100 years (see Fig. 2 of Karmann et al.; 2001). 

2. A more conclusive argument is that the present simula- 
tions show that the presence of ~ 2 IQ^^g of dust in the 
< 10 AU region might be explained, within a moder- 
ate mass disc, by collisional processes alone. Moreover, 
such a low mass disc is consistent with what should be 
expected, in the inner disc, considering the age of the 
system: a disc of debris left over after the accretion 
process of planetary embryos (see next subsection) . In 
short, there is no need for a cometary activity in terms 
of production of the observed dust. 

3. In any case, even if all the dust was to be produced 
by evaporating comets, then the present simulations 
show that mutual collisions within such a ~ 2 10^^ g 
dust disc would anyway be unavoidable. This would 
strongly affect the size distribution, which would prob- 
ably tend towards the steady-state profiles displayed in 
section 4. 

Of course, these arguments are relevant only for the 
inner (3V\c disc. We do not rule out the possibility that 
comet evaporation could be a dominant dust-production 
source in the outer parts of the system, as suggested 



by the Orbital E vaporating Bodies scenario proposed by 
iLecavelieJ ()l998|) for the region beyond 70 AU. 

This might appear to be a somehow paradoxal result, 
since evaporation processes should be more effective in 
the inner regions. But let us once again stress that these 
higher evaporation rates are precisely what makes it diffi- 
cult to find a way to sustain evaporation activity over long 
time scales in the inner disc (as previously mentioned, a 
10 km object evaporates in less than 100 years at 5 AU) . At 
larger distances, volatile evaporation rates are much lower, 
thus reducing the crucial problem of " refilling" the evap- 
oration region with fresh material. Of course, distances 
from the star must remain within the limiting distance 
at which evaporation is possible , i.e. 100-150 AU for CO 
llT.ecavelier et a,].[ ll 99(t. 

5.4. Presence of already formed planetary embryos 

Of course, one cannot rule out the possible presence of 
isolated much more massive objects, such as planets or 
planetary embryos, whose isolation decouples them from 
the collisional cascade resp onsible for the dust pro duction 
(a possibility considered bv lWvatt fc Denti l)2002j) for the 
Fomalhaut system). In fact, the low mass in the dust- 
to-planetesimals range could be interpreted as the conse- 
quence of the presence of such planetary embryos: most of 
the initial mass of the system would already have been ac- 
creted in these embryos, leaving a sparse disc of remnants. 
This would be in accordance with the estimated age of the 
system, a few 10^ years, which significantly exceeds the ex- 
pect ed timespan for the formation of planetary embryos 
(e.g. iLissaueil 1199,31) . so that if embryos have to form, 
then they should be already here. 

Another argument previously proposed in favour of 
the presence of already formed massive embryos is that 
such objects are a good way to explain the disc's thick- 
ness. lArtvmowicd l(l997^ estimated that numerous Moon- 
sized bodies are required in order to induce vertical veloc- 
ity dispersions of smaller bodies of the order o i Q.lvk pr}- 
Nevertheless, as pointed out bv iMouillet et al.l l)l997j) . a 
giant planet on a slightly inclined orbit (like the planet 
required to explain the warp in the outer regions) could 
achieve just the same result: rapid precession of the dust 
particles orbits in the inner regions would lead to thicken 
the disc so that the aspect ratio appears to be equal to 
the planet's inclination. 

6. Conclusions and Perspectives: towards a 
coherent picture of the inner j3 Pictoris disc? 

The present study show that the observed 10^^ to 10^^ g of 
dust in the inner disc is compatible with what would be 
expected from a collisional cascade within a disc of a few 
10-2 M® bodies ranging from micron to kilometre-sized 
objects, without the need for any additional cometary- 
evaporation activity. And even if there was such a 
cometary activity, the required density of objects would 
inevitably lead to important collisional effects. 
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Simulations also show that the size distribution set- 
tles towards a quasi-equilibrium state that strongly de- 
parts from the classical dN oc R~^'^dR Dohnanyi power 
law. This is particularly true for the smallest grains close 
to the radiation pressure ejection limit. However, these 
departures do not too strongly affect the global Dust- 
to-Planetesimal mass ratio and cannot account for the 
incompatibility between the small amount of observed 
dust and the huge number of kilometre-sized FEBs re- 
quested to sustain the transient absorption features activ- 
ity. Furthermore, our runs show that this requested mass 
of FEBs leads to a much too rapidly coUisionaly eroding 
disc that cannot survive on long timescales. The FEB sce- 
nario thus cannot hold in its present form and has to be 
seriously revised. 

We might thus converge towards a coherent picture of 
the inner /3 Pictoris disc: This inner disc should be bound- 
ered by one giant planet of ~ 1 Mjup, located around 
10 AU on a slightly inclined (in order to explain the ob- 
served outer warp) and possibly eccentric orbit (in order to 
trigger the FEB activity). The observed amount of dust 
should be produced by coUisional erosion within a low 
mass disc. Such a low mass disc could be made of debris 
leftover after the accretion of one or several planetary em- 
bryos, the presence of which is fully compatible with the 
estimated age of the system, i.e. a few lO'' years. In other 
words, we should be now witnessing a planetary system in 
a late or at least intermediate stage. The bulk of the ac- 
cretion process is over, but a consequent disc of remnants 
is still present and coUisionaly eroding. 

Our results would also help putting new constraints 
on the SED fits that are usually performed to derive dust 
densities and radial distributions from observed spectra. 
Let us recall that the dust mass estimations for the in- 
ner disc, which we used as inputs for our simulations, 
have been computed e ither by postulating t hat grains 
are of cometary origin ((Li fc Greenberel Il998l) or by do- 
ing a pure mathema tical fit with several free parameters 
i Pantin et al.Lll997(l . In this respect, it would b e interest- 
ing to perform a work similar to that of Li fc Greenberd 
l|l998l) but with a population of coUisionaly produced 
grains as input. An interesting attempt at doing such a 
kind o f study has been recently made bv lWvatt fc Dent I 
in their very detailed study of the Fomalhaut's de- 
bris disc. Nevertheless, their precise fit of the SED was 
made assuming a single power law for the size distribution 
(even though the authors were fully aware of the fact that 
such an academic distribution cannot hold for the small- 
est grains because of the cutoff effect). Such an SED-fit 
analysis goes beyond the scope of the present paper and 
requires additional work. 

It requires in particular to model the whole (3 Pictoris 
disc and not only the innermost parts that only partially 
contribute to the total flux. Only then could the obtained 
size distribution be compared to SEDs integrated over the 
whole disc. A crucial problem would probably be to see 
if the underabundance of millimetre-sized objects that we 
obtained in the inner disc is also to be found for the sys- 



tem as a whole; this would then contradict previous es- 
timates stating that the observed mass of millimetre ob- 
j ects is in accordan ce with a —3.5 equilibrium power law 
l|Artvmowic^ Il997|) . Such a study should of course also 
address more deeply the question of the physical nature 
of the dust grains. It will be the purpose of a forthcoming 
paper. 
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